Time dependent impurity in ultracold fermions: orthogonality catastrophe and beyond 
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Recent experimental realization of strongly imbalanced mixtures of ultracold atoms opens new 
possibilities for studying impurity dynamics in a controlled setting. In this paper, we discuss how the 
techniques of atomic physics can be used to explore new regimes and manifestations of Anderson's 
orthogonality catastrophe (OC), which could not be accessed in solid state systems. 

Specifically, we consider a system of impurity atoms, localized by a strong optical lattice potential, 
immersed in a sea of itinerant Fermi atoms. We point out that the Ramsey interference type 
experiments with the impurity atoms allow one to study OC in the time domain, while radio- 
frequency (RF) spectroscopy probes the OC in the frequency domain. The OC in such systems is 
universal for all times and not only in the long time limit and is determined fully by the scattering 
length and the Fermi wave vector of the itinerant fermions. We calculate the universal Ramsey 
response and RF absorption spectra. In addition to the standard power-law contributions, which 
correspond to the excitation of multiple particle-hole pairs near the Fermi surface, we identify a 
novel important contribution to OC that comes from exciting one extra particle from the bottom 
of the itinerant band. This contribution gives rise to a non-analytic feature in the RF absorption 
spectra, which shows a non-trivial dependence on the scattering length, and evolves into a true 
power-law singularity with universal exponent 1/4 at the unitarity. 

We extend our discussion to spin-echo-type experiments, showing that they probe more compli- 
cated non-equilibirum dynamics of the Fermi gas in processes in which an impurity switches between 
states with different interaction strength several times; such processes play an important role in the 
Kondo problem, but remained out of reach in the solid state systems. We show that, alternatively, 
the OC can be seen in the energy counting statistics of the Fermi gas following a sudden quench 
of the impurity state. The energy distribution function, which can be measured in time-of-flight 
experiments, exhibits characteristic power-law singularities at low energies. 

Finally, systems in which the itinerant fermions have two or more hyperfine states provide an 
even richer playground for studying non-equilibrium impurity physics, allowing one to explore non- 
equilibrium OC and even to simulate quantum transport through nano-structures. This provides a 
previously missing connection between cold atomic systems and mesoscopic quantum transport. 



PACS numbers: 47.70. Nd, 67.85.-d,72.10.-d, 05.60.Gg 
I. INTRODUCTION 



Interest in nonequilibrium quantum dynamics has in- 
creased dramatically in the last few years following exper- 
imental realizations of synthetic many-body systems with 
ensembles of ultracold atoms^ 2 4 With ultracold atoms 
it is not only possible to prepare microscopic systems 
with desired many-body Hamiltonians but, crucially for 
studying dynamics, parameters of such Hamiltonians can 
be changed on time scales that are much faster than 
intrinsic microscopic timescales. Ultracold atomic en- 
sembles are also very well isolated from the environ- 
ment, so states prepared out of equilibrium can undergo 
quantum evolution without relaxation or loss of coher- 
ence^. Finally a rich toolbox of atomic physics makes 
it possible to provide detailed characterization of many- 
body systems, which is crucial for describing complicated 
transient states resulting from non-equilibrium dynamics. 
Recent experimental studies addressed such questions as 
relaxation of high energy metastable states*"^, frydrody- 



namic expansion of strongly interacting fermions in opti- 
cal lattices^, decoherence of split condensates®^, coherent 
su perex change- mediated spin dynamics^, spinor d ynam- 
ic d 11 * 12 !, relaxation and thermalization in ID systens^^, 
as well as interaction quenches in fermionic systems^. 

In this paper we discuss how one can use currently 
available experimental tools of ultracold atoms to study 
the problem of the orthogonality catastrophe (OC) in 
many-body fermionic systems. The core of this problem 
is understanding the response of a Fermi gas to a sud- 
denly introduced localized impurity. This question was 
originally considered in the context of the X-ray absorp- 
tion spectra in metals^, where it manifests itself in the 
characteristic power-law threshold singularities. 

Being one of the very few known examples of a non- 
equilibrium solvable many-body problems, the OC also 
provides a conceptual framework for understanding sev- 
eral fundamental phenomena in solid state physics, in- 
cluding the Kondo effect^', resonant tunneling in meso- 
scopic structures" 2 ^ 1 , ID quantum physics beyond the 
Luttinger liquid paradignP3H22 an( j the motion of a 
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FIG. 1. (Color online) Proposed experimental setup: Impu- 
rity atoms (blue dots with arrows, indicating internal states) 
are immersed in a sea of itinerant host fermions (red dots). 
The optical lattice, indicated by the black parabolas, localizes 
the impurity atoms without affecting itinerant host fermions. 
We assume a low concentration of impurity atoms, so that we 
can consider scattering on a single impurity. The scattering 
of fermions can be controlled by applying magnetic fields and 
by manipulating the internal impurity states with RF-pulses. 
The system's parameters can be controlled quickly (compared 
to the intrinsic timescales of the many-body system of host 
fermions), which makes this setup ideal for studying non- 
equilibrium impurity physics, including the new regimes of 
OC. 



heavy particle in a Fermi gas^. 

The basic setup that we focus on is shown in Fig.[T] 
We consider a quantum degenerate mixture of two types 
of atoms, e.g. a Bose-Fermi mixture^H2H or a Fermi- 
Fermi mixture^HUl. We assume that one type of atoms, 
called an impurity atom below, has a much lower den- 
sity than the other. The majority atoms, referred to as 
host atoms, are always taken as fermionic. Two types 
of atoms can have very different polarizability hence it 
is possible to create an optical lattice that strongly lo- 
calizes impurity atoms while having very little effect on 
the host fermions. The hyperfme spin states |t) , \i) of 
the impurity atoms can be manipulated using RF fields, 
which allows to switch between weakly and strongly in- 
teracting states with respect to host fermions and thus 
introduce time dependent impurities in the Fermi ga:P^. 
Additionally, the strength of interactions between host 
atoms and a given hyperfme state of the impurities can 
be controlled by tuning magnetic fielcP^Hifl. We assume 
a regime of very low density of impurity atoms so that 
scattering processes taking place on different impurities 
can be analyzed separately. Thus in the rest of the paper 
we will discuss dynamics of a single impurity interacting 
with a Fermi gas3S. 

The setup proposed above, with its tunability and pre- 
cise control, provides a way to study new regimes and 
manifestations of OC, which remained out of reach in 
solid state systems. First, we will show that the OC is 



a fully universal function of the Fermi wave vector and 
the impurity-host-fermion scattering length; this func- 
tion will be calculated below. The full solution reveals 
new singularities in the RF absorption spectra, which 
emerge away from the absorption threshold. Second, we 
will demonstrate that the Ramsey and spin-echo spec- 
troscopy - standard tools used in cold atoms - pro- 
vide a direct access to the OC, as well as to more com- 
plex non-equilibrium response functions, in the time do- 
main. Third, manifestations of OC in the energy count- 
ing statistics will be discussed. Finally, by generalizing 
the setup described above, it is possible to study OC in 
the regime when the Fermi gas is driven out of equilib- 
rium, as well as to simulate quantum transport. 

Given the lon g illu strious history of the original prob- 
lem, see Ref s ESHsiHSH&il it is useful to start by summa- 
rizing the new aspects of OC in cold atomic systems. 

Extended universality. Experiments with ultracold 
atoms can demonstrate the universality of OC in a much 
broader sense than it was previously discussed for elec- 
tron systems. The most basic quantity demonstrating 
OC is the time-dependent overlap function, 

S(t) = (ikle^We-* 6 **'*^), (1) 

where \ipo) is the initial ground state of the Fermi gas, 
and Hi<t) is the Hamiltonian before(after) the impurity 
is switched on. In solid state systems the universality 
was understood only as a statement about the long time 
behavior of S(t). At long times, t ^> h/Ep, Ep being the 
Fermi energy, the overla p ex hibits a power-law decay^l, 
reflecting Anderson's OC 51 * 52 l The power-law decay is 
the result of the large phase space for exciting multiple 
low-energy particle-hole pairs following the switching of 
the impurity. 

Although OC has been considered as an "exactly solv- 
able" dynamical many-body problem, in solid state sys- 
tems model parameters are either not known accurately 
or provide a crude approximation to much more com- 
plicated many-body processes. For example, a local im- 
purity potential comes from the Coulomb potential of a 
hole screened by the conduction electrons. A one particle 
scattering picture commonly used in the analysis of the 
OCPHSU, which neglects many-body aspects of electron- 
electron interactions, is only an approximation valid for 
low energy scattering processes of electrons close to the 
Fermi energy. Furthermore, the density of states for the 
conduction band electrons is typically not known and can 
be modified by the electron interactions. In addition, ex- 
trinsic degrees of freedom, including phonons, can play 
an important role in the OC dynamics leading to extra 
decay factors. Hence one can neither calculate the full 
time dependence of S(t) nor claim its universality at all 
timescales. 

By contrast, in the case of ultracold atoms, one cre- 
ates a gas of Fermi atoms, which are well isolated, do not 
interact with each other, and have a simple parabolic dis- 
persion. The interaction of the localized impurity atoms 
and conduction band fermions is fully characterized by 
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a single parameter, the scattering length a. In the case 
of wide resonances, one finds a universal behavior of the 
energy dependence of the scattering amplitude 



1 



1 .V2mE 



(2) 



where E is the energy of the scattering atoirPH. Such 
scattering amplitude leads to the following universal be- 
havior of the scattering phase shift on the momentum of 
the scattered fermion: 



8(k) — —tan x ka, 



(3) 



where — tt/2 < S(k) < tt/2. Consequently, the entire 
function S(t) (and not only its long time asymptote) is a 
universal function of kpa and Ept (hp is the Fermi wave 
vector and Ep is the Fermi energy), 



S(t) = S(k F a,E F t/h). 



(4) 



The main result of this work is the calculation of universal 
functions S(t), as well as more complicated time-domain 
response functions. Specific examples are shown in Figs. 
[3j§and[7] 

Measuring orthogonality catastrophe in time domain. 
In solid state systems the OC ([!]) is observed indirectly 
in the frequency dependence of the absorption spectrum, 
given by 



1 f c 
A u = -Re / 

T Jo 



f 5(t) dt. 



(5) 



The power-law decay of the overlap function in the time 
domain translates into p ower-law threshold singularities 
in the absorption spectra^- 6 * 50 * 52 !. In ultracold atoms one 
can perform similar measurements of OC in the frequency 
domain by doing RF spectroscopy on impurity atoms. 
Examples of the universal absorption spectra for different 
values of kpa are shown in Fig.[4j 

However, it may be more illuminating to measure S(t) 
in time domain using Ramsey type interference, which 
is a well established experimental technique in atomic 
physics. While it was originally designed for metrology 
applications, it has been realized recently t hat it can also 
be used as a probe of many-body dynamic d 8 * 57 * ^. 

The idea of the proposed experiment is as follows: ini- 
tially, the impurity is prepared in the down-state |J,), and 
the fermions are in the ground state \ipo}- Then, the 
Ramsey interferometry is performed: at time t = 0, a 
tt/2 pulse is applied, such that the system is driven into 
the superposition state <S> \ijja). The two states 

evolve differently since the |f) and \l) states interact dif- 
ferently with the Fermi sea: 



!*(*)) = \i) ® e-^/ R |Vo) + It) ® e-^/^o). 

The fermions stay undisturbed in the first state, while the 
impurity scattering excites multiple particle-hole pairs in 



the second state. Performing a second tt/2 pulse after 
time t, and measuring (^>(t)\S x \^f(t)} gives: 



(Sx(*)> 



Re{iPo\e iHit/n e- iH f t/n \ip ) 



ReS(t). (7) 



In the equation above we neglected the trivial phase fac- 
tor arising from the energy difference of states |f) and 
I4-). Thus, the Ramsey interferometry provides a direct 
measurement of the OC overlap^. 

The basic Ramsey type experiment described above 
corresponds to a local quench type dynamics in which 
the impurity strength is changed once. One can also use 
the Hahn spin echo, as well as more complicated spin- 
echo-type sequences, to realize processes in which the 
impurity strength effectively switches between different 
values multiple times. As we show below the response 
of the Fermi gas to such processes is characterized by 
a power-law decay of the overlap, with an exponent en- 
hanced compared to the case of usual OC. 

The predicted faster decay of the spin-echo response 
goes against the atomic physics intuition; it stems from 
the fact that the spin-echo sequence does not "undo" the 
evolution of the Fermi gas under impurity scattering, but 
instead, creates a state in which the impurity pseudospin 
and Fermi gas are strongly entanglcd 63 4 From the exper- 
imental point of view spin echo experiments should have 
an additional advantage that they allow to cancel slow 
fluctuations of the magnetic field. 

New universal features: bottom of the band physics. 
What are the new universal characteristics that one can 
find in S(t), as well as in the corresponding absorption 
spectra A(u) at intermediate time scales? One feature 
that can be seen from Fig. [3] is the oscillations of S(t) on 
the timescale HEp 1 . The origin of this novel feature is 
the following: after the impurity is introduced, there is 
a class of processes in which, in addition to a large num- 
ber or low-energy particle-hole pairs, an extra hole with 
energy ~ Ep is excited near the bottom of the band. 
The phase space of such processes is enhanced, owing 
to the van Hove singularity in the density of states at 
the band bottom in one dimension (note that OC is es- 
sentially a ID problem since only the s-wave channel is 
important). The combined dynamics of the high-energy 
hole and the low-energy particle-hole pairs result in an 
additional weaker power-law contribution to the overlap 
function. 

This contribution to the overlap function becomes even 
more evident in the frequency domain, where it gives rise 
to a cusp singularity at the energy Ep above the thresh- 
old (see Fig. [3]). We find that, as the Feshbach resonance 
is approached from the side of the negative scattering 
length, the cusp develops into a true singularity with a 
universal exponent 1/4; for any finite value of kpa, this 
peak is smeared, but its overall shape is described by a 
universal function which is discussed below. We note that 
such behavior of the response function across the Fesh- 
bach resonance is rather unusual, since it shows a truly 
singular behavior only exactly at the resonance. For com- 
parison, in conventional BCS-BEC crossover studiesP, 
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most of the measured quantities show a smooth behavior 
exactly at the resonance. 

Non-equilibrium orthogonality catastrophe. So far we 
assumed that the itinerant host fermions do not have any 
internal degrees of freedom. Generalization of the setup 
shown in Fig. [T] to the case of a multi-component Fermi 
gas (with multiple internal states) allow one to realize 
a wide range of dynamical impurity phenomena in non- 
equilibrium Fermi gases. This is an even richer class of 
problems that arises mesoscopics^ in particular, prob- 
lems of this kind describe the quantum transport through 
any mesoscopic structure (e.g., a point contact), where 
the Fermi seas in two or more leads are kept at different 
chemical potentials. Mathematically, such problems can 
be reduced to a non-abelian Riemann-Hilbert problenPSl, 
the solution of which is not known in a general case. 

While usually transport is difficult to study in systems 
of ultracold atoms (see however^), we will show that 
the extension of our setup provides a way of simulating 
the quantum transport. More generally, it allows one to 
study the response of the non-equilibrium Fermi gas in a 
controlled setting. 

New quantum observables and distribution functions. 
OC experiments with ultracold atoms should make it pos- 
sible to study quantum variables that are not accessible 
in electron systems. For example, the full ener gy of an 
interacting many-body system can be measure d 46 * 67 * ^. 
Moreover, it is possible to measure not only the average 
values but also fluctuations of quantum observa bles-^, 
and in some cases even the full distribution function d 61 l 71 l 

In the quantum impurity system as in Fig.[T] following 
the impurity potential quench the system is no longer in 
an energy eigenstate and the full distribution function 
of the total energy should also exhibit power-law type 
singularities (see Refs. [201 and 172] and discussion below), 
which provides an alternative manifestation of the OC. 

Full counting statistics of scattering processes should 
provide an intriguing connection to an extensive theoret- 
ical research in mesoscopic physics^. It is worth noting 
that the measurements of charge counting statistics were 
notoriously difficult in solid state systems^. Given the 
available experimental tools, such measurements should 
become possible in cold atomic systems, with an addi- 
tional advantage that the full counting statistics of scat- 
tering events can be measured for fermions in specific 
energy windows rather than in the whole energy range. 

This paper is organized as follows: In Sec. [Tl] we de- 
scribe our approach to OC, and present the results for 
the universal overlap function S(t), which can be mea- 
sured by the Ramsey interferometry. The universal radio- 
frequency (RF) spectra A{uS) are evaluated, and their 
properties, singularities, as well as the novel "bottom of 
the band" feature are discussed. Extensions to the Hahn 
spin-echo and more complicated spin-echo type experi- 
ments, in which effectively the impurity strength changes 
several times, are studied in Sec. |III| The manifesta- 
tions of the OC in the energy counting statistics and 
generalizations to the nonequilibrium OC using multi- 
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FIG. 2. (Color online) Schematic illustration of the inter- 
mediate excited sets of states giving power-law contributions 
to the overlap function, which result in the threshold sin- 
gularities in the RF absorption spectrum, (a) Undisturbed 
Fermi sea (impurity is in the state), (b-c) Case a < 
(no bound state): (b) Multiple low-energy particle-hole pairs, 
giving rise to the leading power-law decay ( |12[ ) of the overlap; 
(c) an additional particle is promoted from the bottom of the 
band to the vicinity of the Fermi surface. This contribution 
leads to the oscillations with period 2hTi/Ef in the overlap 
function (see Fig.|3^a)), as well as to the new cusp-like singu- 
larity in the absorption spectra (Fig.|4|a)). (d-f) For the case 
a > 0, when the impurity potential creates a bound state, 
there are three important sets of states: (d) the bound state 
is empty and multiple low-energy excitations are created, (e) 
the bound state is filled, and (f) an additional particle is ex- 
cited from the bottom of the band to the Fermi surface. The 



processes (d) and (e) lead to the behavior ( 15 1 of the overlap 
function, while (f) leads to a faster decaying oscillating con- 
tribution. The three contributions result in the absorption 
spectra shown in Fig. |4jT). 



component host atoms are discussed in Sec. IV and Sec. 
|Vj respectively. Finally, we connect the proposed setup 
to existing cold-atom experiments and conclude our find- 
ings in Sec. rVT] 



II. 



UNIVERSAL ORTHOGONALITY 
CATASTROPHE 



In this Section, our main goal will be to calculate the 
universal response functions, S{t) and A{u). To sim- 
plify the discussion we assume that the \l) state of the 
impurity does not interact with host fermions and only 
the |t) state gives rise to scattering. Generalization of 
our analysis to the case when both |f) and ||) states in- 
teract with the fermions is straightforward. Assuming 
that the impurity atom is in the rotationally symmet- 
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ric ground state of the confining parabolic potential and 
cannot be excited to higher states (the energy of the con- 
fining potential is much larger than the Fermi energy), 
we only need to consider the s-wave scattering of host 
atoms on the impurity. The corresponding Hamiltonian 
can be written in the following form: 



H = H + |t) (tl ® V, 



(8) 



where Hq = SkC k Ck is the Hamiltonian of free 

fermions, and V = Vq J^k q c l c <? describes scattering. 
Physically, the scattering potential is characterized by 
the scattering length a. The Ramsey sequence described 
above measures the real part of the overlap function, see 
Eq.(B, where Hi = H Q and H f =H + V. 

There are two physically distinct cases, for which the 
asymptotic behavior of the overlap function (as well as 
the shape of the absorption spectrum) is qualitatively 
differenpS (i) a < 0, when the impurity potential does 
not create a bound state, and (ii) a > 0, when there is a 
bound state with energy 



E h = 



2ma 2 



(9) 



To gain an intuition about the behavior of the overlap 
function S(t) in the two cases, it is convenient to rewrite 
it by inserting a complete set of eigenstates {ip a } (with 
corresponding eigenenergies {E a }) of the Hamiltonian 
H + V into Eq.|l]): 



2 e i(E -E a 



)t/h 



(10) 



In Fig. [2] we illustrate the dominant intermediate states 
\ip a ) giving rise to the powerlaw contributions to the over- 
lap. In the case a < 0, the main contribution to S(t) is 
due to the intermediate states \ip a ) in which multiple 
low-energy particle-hole pairs are created near the Fermi 
surface (see Fig. |2jb)). The large phase space of such 
excitations gives rise to the power-law decay of S(t) at 
long times, see Eq.([l2|). In the case a > 0, when there 
is a bound state, there are two distinct important sets 
of states \ip a )i which involve many excitations near the 
Fermi surface, but differ by the bound state being ei- 
ther empty or filled (Fig. |2jd),(e)). The contributions of 
these two sets of states into S(t) separate, such that at 
long times S(t) is given by the sum of two power-laws 
with different exponents, see Eq.(15). 



Below we identify another important set of interme- 
diate states, which gives rise to an additional weaker 
power-law contribution to overlap. Such states involve 
a single hole excited from the vicinity of the band bot- 
tom, in addition to a number of low-energy particle-hole 
pairs (Fig. [2^c,f ) ) ; the contribution of these states to the 
overlap function is enhanced by the van Hove singularity 
at the bottom of ID band. As we shall see below, this 
new contribution results in an unusual kink feature in the 
absorption spectrum, which at the unitarity (|/cfq| = 00) 
develops into a full power-law singularity with a universal 
exponent 1/4. 




FIG. 3. (Color online) Universal overlap function |5*(i)|, 
which can be measured in the Ramsey interferometry experi- 
ment, is shown for different values of the scattering length a, 
see legend. Panel (a) shows \S(t) \ for the case a < 0, when the 
impurity potential does not create a bound state. The dashed 
lines are power-law fits ( 12 1 at long times. At intermediate 



times, oscillations coming from the bottom of the band are 
visible. Panel (b) shows \S(t)\ for the case a > 0, when the 
impurity potential creates a bo und state. The asymptotic 
behavior is described by Eq.(15l, and exhibits strong oscilla- 
tions. 



A. Method 

Our analysis of the universal OC is based on the repre- 
sentation of the respon se functions in terms of functional 
determinant d 20 * 21 * 53 * 74 !. We are interested in evaluating 
many-body averages of products of exponents of oper- 
ators quadratic in creation/annihilation operators [see 
Eqs. ^ and Q]. Such many-body quantities can be re- 
duced to functional determinants in the space of single- 



G 



particle orbitals 
written in the following for 



The overlap S(t), 



ap 



20I21I74L 



in particular, can be 



S(t) = det [1 - h + R(t)h) , R(t) = e iho*/* e -*V*o+1>)t/h > 

(11) 

where h is the occupation number operator, n\e) = 
n(e)|e); at finite temperature, n(s) is given by the 
Fermi distribution, n(e) = exp (( e _)j)/T) + i ■ The opera- 
tors ho, ho + v are the single-particle Hamiltonians in 
the absence and presence of im purity, respectively. 

While previous workPSK^SJZSI concentrated on analyz- 
ing the asymptotic behavior (at times t ^> h/Ep) of the 
functional determinant of the type (111, here our goal is 



to find the response functions at all times. Therefore, we 
numerically evaluate the functional determinants. We 
consider the case of a finite system confined in a ball. 
This allows us to view the operator in Eq.(ll| as a 



finite-dimensional matrix for the case of zero tempera- 
ture, when n(e) is non-zero only for a finite number of 
states under the Fermi level. It turns out that the ma- 
trix elements of the matrix R(t) can be evaluated to high 
precision by appropriately truncating the (infinite dimen- 
sional) matrix e - l ( h o+v)t _ Taking the finite size scaling, 
we obtain the universal behavior of the overlap in the 
continuum limit. The details of our discretization and 
truncation procedures can be found in Appendix [X] An 
advantage of our approach is that it can be easily gener- 
alized to other geometries, including the experimentally 
important case of the harmonic trap. 

The structure of the single-particle Hilbert space de- 
pends on the sign of the scattering length: when a < 0, 
the eigenstates of the Hamiltonian ho + v are scatter- 
ing states; when a > 0, the spectrum should be supple- 
mented by the bound state. This translates into a differ- 
ent asymptotic behavior of the overlap function S(t) and 
the absorption function A(u)), as discussed below. 



B. Universal overlap functions 

The absolute value of the resulting universal function 
S(t) in the two cases is shown in Fig.[3j (We note that the 
quantity Re S(t), measured in the Ramsey experiment, is 
quickly oscillating with a period set by the energy differ- 
ence between two hyperfine states as well as by the dif- 
ference of the ground state energies of the Hamiltonians 
Hq and Hq + V; when plotting the Ramsey response in 
Fig. |3j we have ignored this trivial phase and illustrated 
the Ramsey amplitude |iS'(t)|-) 

In the case a < 0, see Fig.[3](a), the long-time asymp- 
totic behavior at t ^> h/Ep agrees with the analytic re- 
su! 



S(t) « Ce- lAEt/h 



1 



where 



iEpt/h + 



-tan 1 (kpa) 



(12) 



(13) 



k a=-0.5 



k p a=-1 .5 
k p a=-3.0 
k p a=-6.0 




FIG. 4. (Color online) RF absorption spectra for different 
values of scattering length a. (a) When a < no bound 
state is present and the spectrum exhibits one power-law 
edge singularity at TkJo = AE, due to the processes illus- 
trated in Fig. [2^b). In addition, a weak non-analyticity at 
huii — AE + Ef is found, attributed to the excitations from 
the bottom of the band, illustrated in Fig. [2^c). This non- 
analyticity is magnified and plotted for additional large val- 
ues of fcpa = — 12, —18, —24 in the inset in (a). At unitarity, 
the non-analyticity transforms into a true divergence with ex- 
ponent 1/4, described by a universal result (181. (b) When 
a > 0, the impurity creates a bound state. The spectra in 
this case exhibit two threshold singularities, which are offset 
by an energy \Ef — Eb\, as well as an additional cusp-like 
singularity the origin of which is illustrated in Fig. [2jf ) . 



is the scattering phase shift at the Fermi wavevector, 



AE 



dE 



S(V2mE) 



(14) 



is the energy renormalization of the Fermi sea due to 
the impurity leveP^. The power law in Eq. ( 12 ) is a 
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manifestation of the "shake-up" process after the sudden 
switching of the impurity potential, which creates mul- 
tiple low-energy particle-hole excitations near the Fermi 
surface, see Fig. 2jb). The dependence of the prefactor 
C on kpa, calculated from fitting the numerical results 
at long times with formula (12), is illustrated in Fig. [5] 
The overlap function S(t) exhibits oscillations with pe- 
riod 2nh/Ep (see Fig. [3]), which are due to the processes 
in which an extra particle is excited from the bottom of 
the band to the Fermi level (see Fig. [2^c) ) . 

For the case a > 0, see Fig.[3](b), the overlap func- 
tion S(t) has two main contributions, coming from the 
intermediate states where the bound state is cither filled 
or empty. This alters the asymptotic behavior, which is 
now represented by the sum of two power laws^ 



S{t) ^Ce- tAEt/h 



1 



iE F t/h + 



q -i{AE-E F +E h )t/H 

a b = (l + 5 F /ir) 2 , 



iE F t/h + 



(15) 



where C and Cb are fcj?o dependent numerical constants. 
The strong oscillations with period 2ttTi/(Ep — Eb), pre- 
dicted by the above formula, are evident in Fig.[3jb). Our 
approach allows the numerical calculation of the prefac- 
tors C and Cb, yielding the result shown in Fig. [Bp 
In addition, the overlap function exhibits faster-decaying 
oscillations with period 2Trh/Ep, which correspond to 
the "bottom of the band" contribution (see Fig[2]^f)); in 
Fig. [3jb) , these are not visible because they are masked 
by the stronger oscillations with period 2irfr/ '(Ep — Eb). 



C. Universal radio- frequency spectra 

We now use the above results for the overlap func- 
tion to evaluate the RF spectra. Numerically calculating 
S(t) for a very large interval of t values, necessary for 
a precise Fourier transform, is computationally challeng- 
ing. We circumvent this difficulty as follows: we choose 
a cut-off f* 3> h/Ep; at times t < t*, function S(t) is 
evaluated numerically, while at t > i„, we match S(t) to 
its asymptotic power-law behavior given by Eq. ( 12 ) for 



a < and Eq. ( 15 ) for a > 0. Calculating the real part of 



the Fourier transform of the resulting function, we obtain 
the RF spectra. 

The behavior of A(u>) for negative and positive scat- 
tering length is qualitatively different (see Fig. [4j. The 
absorption spectrum for a < exhibits one edge singu- 
larity, with an exponent 1 — a that follows from Eq.( 12 ). 



For a > 0, the asymptotic behavior ( 15 ) gives rise to a 



double-threshold absorption spectrum; due to the pres- 
ence of the bound state the spectrum is characterized 
by two singularities at energies hub = AE — Ep + Eb 
and Tluiq = AE. The first threshold corresponds to fill- 
ing the bound state following absorption, and the second 
to leaving it empty^H As follows from the definition of 
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FIG. 5. (Color online) The prefactors that appear in 
the asymptotic behavior of the overlap function, as well 
as in the RF response near the thresholds at cjo,Wb, see 
Eqs. p2||,|[T5|,|[l6||, and |T7| . The prefactors were obtained 
from the asymptotic form of the numerically evaluated over- 
lap functions. 



A(ui), Eq. |5]), and from the asymptotic form (15), the 
two threshold singularities are characterized by different 
exponents, 1 — ab and 1 — a, respectively. The excitation 
processes which correspond to these singularities are il- 
lustrated in Figs. [2je) and[2]jd). The respo nse nea r these 
thresholds can be obtained using formulas ( 12|15 1, which 
gives: 



S(u- loq) 
S(oj- u b ) 



2irC9(u> - lu ) 
r(a)(E F /h)<* 



|w-wo| a_1 , 
hu = AE, (16) 

huj b = AE-Ep + Eb. (17) 



The most interesting new feature of the universal spec- 
tra shown in Fig. [4] is the non-analyticity of A(oj) at fre- 
quency huJi = AE + Ep for a < and huji = AE + E b 
for a > 0; the physical origin of this feature was al- 
ready discussed above. As the unitarity is approached, 
the non- analytic feature becomes more prominent. This 
phenomenon, as well as the full structure of A{uS) at 
fcp|a| 3> 1 can be understood as the result of a non-trivial 
interplay between two-body physics that involves the im- 
purity and the hole near the bottom of the band, and the 
dynamics of multiple particle-hole excitations created at 
the Fermi surface. The possibility of a non-trivial in- 
terplay between many-body and few-body physics is a 
unique feature of ultracold atom physics, and has at- 
tracted significant theoretica l 78 * 79 ! and experimental^ in- 
terest recently. 



We have developed an analytic theory of this new fea- 
ture; here, we briefly summarize the essence of our ap- 
proach and main results, providing the full results in Ap- 
pendix [B] the details of the solution will be provided in 
Ref. [SI] - The idea is that, in the time domain, the contri- 
butions from the Fermi surface excitations and from the 
single hole at the bottom of the band to the determinant 
(111 factorize. The former contribution is given by the 
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power law which can be obtained within the functional 
determinant approach, while the latter can be evaluated 
exactly within the two-body theory^. 

In the frequency domain, the result can be written as 
a convolution of the two terms describing these two pro- 
cesses (see Appendix [B]); this gives rise to a cusp-like sin- 
gularity at an energy AE + Ep for a < and AE + Ei, 
for a > 0. Away from the unitarity, this singularity is 
smeared on the energy scale ~ h 2 /2ma 2 (for either sign 
of the scattering length a); the origin of this smearing 
lies in the dynamics of the hole. 

Remarkably, right at the unitarity, \kpa\ — oo, 
S 2 /2mo 2 vanishes, and a full non-smeared power-law sin- 
gularity with the universal exponent 1/4 develops; this 
singularity is asymmetric, and is described by: 



A(co) 



1.74^-cj!!- 1 / 4 
(£fA) 3/4 

r(i/2) 



(cj - UJi) 



r(3/4) 



+ 0(wi - lu) 



r(i/4) 
r(i/2) 



(18) 



for \kpa\ = oo and h\u> — uj\ \ <C Ep. 

When a < 0, this peak gets smeared out at energies of 
the order of h 2 /2ma 2 from its maximal value; the evolu- 
tion of the peak depending on the value of kpa is illus- 
trated in the inset to Fig.[4ja). 

For the case a > 0, the true bound state with energy Ej, 
"pinches off' from the bottom of the band and leads to a 
threshold with an exponent 3/4. Thus for large but finite 
kpa > the universal form of A(uj) near lo\ has a char- 
acteristic double peak structure, as is seen in Fig. [5|b) 
for kpa = 6.0. For increasing interaction parameter kpa 
the bottom of the band feature approaches the singular- 
ity at wo where the bound state is empty, as the energy 
difference is only of the order of Ef, which tends to zero 
when the unitarity is approached. 

The universal contribution due to the bottom of the 
band is a unique feature of ultracold spinless fermions 
which has no analog in conventional solid state systems. 
Indeed, the universal behavior coming from the excita- 
tions in the vicinity of the Fermi surfaces is ubiquitous 
in solid state systems, and is controlled by Fermi liq- 
uid theory^. On the other hand, away from the Fermi 
surface fermionic quasiparticles are not well defined in 
3D interacting systems, and thus bottom of the band 
contributions to the orthogonality catastrophe cannot be 
probed in solid state systems. In contrast, for spinless 
cold atoms, the bottom of the band contributions are 
well defined, since the interactions between atoms in the 
s-wave channel are absent. Thus, fermionic excitations 
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FIG. 6. (Color online) The spin-echo response ( 19 1 of the 
Fermi gas. At long times, it is characterized by a power-law 
decay ( 20 1 with an exponent three times larger than that of 



the standard OC. 



are well defined for all energies, including the vicinity of 
the band bottom. 



III. SPIN ECHO RESPONSE 

Now we proceed with discussing spin-echo-type se- 
quences, which allow one to measure the response of the 
Fermi gas to a process in which the internal impurity 
states switch multiple times. The spin-echo sequences 
have an important advantage compared to the Ramsey 
intcrfcromctry in that they are not sensitive to the slowly 
fluctuating magnetic fields typically present in experi- 
ments. 

The main effect of the fluctuating fields is to introduce 
an energy difference As between the |f) and ||) states. 
Generally, this would affect the outcome of the Ramsey 
experiment: owing to the different phase picked up by 
the |t) and ||) states, the measured value of (S x ) is given 
by 



S(t)e 



iAst/(2h) 



S*(t)e 



-iAet/(2h) 



<&(*)) = 



The phase factor e lA6t /( 2h ) [ s different from shot to shot; 
thus, averaging over different shots gives rise to an addi- 
tional decay of (S x (t)). This extra decay would compli- 
cate the observation of the OC. Similarly to quantum op- 
tics experiments, the undesired contribution of the fluc- 
tuating magnetic fields can be eliminated in the spin-echo 
experiment. 

We first consider the Hahn spin-echo protocol: ini- 
tially, the system is prepared in the state ||) <g> |^o)j at 
time —t a ir/2 pulse is applied, followed by the ir pulse 
at t = 0. Finally, at time t another tt/2 pulse is applied. 
Similar to the case of the Ramsey experiment, we can 
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FIG. 7. (Color online) The generalized spin-echo response (21 1 of the Fermi gas for n — 1,2,3 compared to the Ramsey 



interference. At long times, the generalized spin-echo response is characterized by a power-law decay (22 1 with an exponent 
(4n — 1) times larger than that of the standard OC. The gray lines correspond to T = 0, while the colored lines are for finite 
temperature (a) T = 0.01 Tf, (b) T = 0.03 Tf, and (c) T = 0.05 Tf. The scattering length is k F a = -1.5 in all cases. The 
power-law decay of the spin-echo response gives way to a faster exponential decay at times t > Ti/T. 



show that S x measured after such a sequence is given by: 
(S x , 1 (t))=Re[S 1 (t)],S 1 (t) = (^ol^Wo), (19) 
where 

(j x _ e -iH t/H e -i(H +V)t/H fj^ _ e -i(H +V)t/h e ~iH t/h 

are the operators describing the evolution of the Fermi 
gas state that was initially coupled to the |f) and |4_) 
states of the impurity. This correlation function describes 
non-trivial Fermi gas dynamics in a process where the 
impurity switches between different states several times. 
Similar response functions arise in the analysis of the 
Kondo problem 17 . To understand the behavior of the 
spin-echo response, we represent it as a functional deter- 
minant given by Eq.(ll) with R = u^iii (ui,tt2 being 



the single-particle analogues of operators Ui, 1/2). As can 
be shown analytically using the method of Ref. [30] , the 
asymptotic behavior of such a determinant at long time 
t fi/Ep is given by a power-law, 



(20) 



with an exponent three times larger than that of OC. The 
universal behavior of the spin-echo response for arbitrary 
times, calculated numerically using the determinant ap- 
proach, is compared to the standard OC probed with 
Ramsey interference in Fig. [6] 

The predicted faster decay of the spin-echo response is 
somewhat unexpected: the intuition from quantum op- 
tics would suggest that the spin-echo response generally 
eliminates the broadening due to slowly fluctuating en- 
vironment, and therefore should be characterized by the 
slower decay compared to the Ramsey interference. This 
intuition, however, relies on the assumption that the en- 
vironment only affects the relative phase of the two hy- 
perfine states; thus, the operators describing the effect 
of the environment on the two pseudospin states com- 
mute (and therefore the spin echo sequence can "undo" 



the effect of the environment). The above discussion il- 
lustrates that this assumption breaks down for the case 
when the Fermi gas plays the role of environment: the 
reason is that the operators e - % ( H o+V)t an( j e ~iH t nQ 
longer commute, and therefore U\ 7^ U2', this results in 



the non-trivial form of the spin-echo response ( 19 ) and 
its faster decay. 

The generalized spin-echo sequence, in which (2n — 1) 
7T pulses are applied at equal time intervals t (the Hahn 
spin-echo corresponds to n = 1), allows measurements 
of even more interesting response functions. The corre- 
sponding overlap function is given by: 

(&,»(*)) = Re[S n (t)], S n (t) = (^ |l> 2 - n t/f |Vo). (21) 

The long-time asymptotic behavior of the nth response, 
calculated similar to the case of the Hahn spin-echo, is 
characterized by the power-law decay 



S n (t) cx t 



An - 1. 



(22) 



Thus, by increasing the number of pulses n, the expo- 
nent can be strongly enhanced. This should facilitate 
the observation of the power-law decay in the spin-echo 
experiments. The universal spin-echo responses for dif- 
ferent values of n, calculated numerically, are illustrated 
in Fig. [7] along with Ramsey interference results. The 
asymptotic form at large times agrees with the analytic 
formula p2] ). 

Experiments are always done at low, but finite temper- 
ature; therefore, it is important to understand the effect 
of non-zero temperature T on the spin-echo response. 
Within our approach, the finite temperature response 
is found by calculating the corresponding determinants 
with the distribution functions n(e) = cxp ^ £ „j J - ) / r - )+1 ■ 
The result, illustrated in Fig. [7J shows that the finite 
temperature does not affect the power-law behavior of 
the response at relatively short times t < fi/T; however, 
at longer times, t > h/T, the responses are character- 
ized by a faster, exponential decay. It is evident from 
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Fig. [7] that two most interesting features of the spin-echo 
response - the power-law decay and the oscillations with 
period 2irfi/Ep - survive up to experimentally accessible 
temperature T = Q.05Ep, and therefore they should be 
observable with current experimental means. 



The power-law singularity ( 24 1 in the energy fluctuations, 
similar to that found in the context of spin systems by 
SilvcP^l, provides a new experimentally accessible mani- 
festation of OC. 



IV. MANIFESTATIONS OF THE OC IN 
ENERGY COUNTING STATISTICS 

So far, we have considered two ways of observing OC — 
the RF spectroscopy (energy domain), and the Ramsey, 
as well as spin echo sequences (time domain) . Both meth- 
ods are based on studying the impurity properties. How- 
ever, cold atomic system also allow direct measurements 
of the Fermi gas properties (e.g., measuring the occu- 
pation numbers of states with different momenta in the 
time-of- flight experiments). Can one see traces of OC by 
measuring properties of the Fermi gas following sudden 
change of an impurity potential? 

An obvious candidate quantity is the time-dependent 
density profile of the Fermi gas. We have studied the 
time-dependent density profile p(t, r), finding that it does 
not carry signatures of the OC. This is because the over- 
lap function involves a complicated sum of the n-particle- 
hole pair correlation functions, while the density operator 
can create at most one particle-hole pair. Thus, we are 
forced to consider the fluctuations rather than averages 
of physical quantities. 

We note that the OC is formally related to the dis- 
tribution function of the fluctuations of the total energy 
of the Fermi gas following an impurity quench. For sim- 
plicity, let us consider a situation where the impurity is 
initially in the interacting state such that fermions 
are in the appropriate ground state \ipo); at some time, 
the impurity is suddenly flipped to the non-interacting 
state |4,). Then, the characteristic function of the total 
energy distribution (with new Hamiltonian) is given by: 

x(A) = (^|e iA ^|^). 

Comparing with Eq.Q, we see that up to a phase factor, 
the characteristic function is identical to the overlap func- 
tion, with parameter A playing the role of time. Thus, 
x(A) exhibits power-law behavior at large A: 



X (A) cx A- 5 -/* 



(23) 



The power-law asymptotic behavior ( 23 ) gives rise to 



the universal power-law behavior of the probability dis- 
tribution P{E) at low energies E -C Ep\ 



P{E) = 



1 

2^ 



— i\E 



(24) 



The unusual distribution function can be directly mea- 
sured in the time-of-flight experiments (the total energy 
is given simply by the kinetic energy since in the final 
state the fermions are not interacting with the impurity) . 



NON-EQUILIBRIUM OC AND QUANTUM 
TRANSPORT 



The orthogonality catastrophe is modified qualita- 
tivel y wh en the Fermi gas is driven out of equilib- 
riurr l 2 ° l 21 4 In solid state systems, such a situation oc- 
curs naturally when an impurity is coupled to two ID 
leads with different chemical potentials. The impurity 
state controls the transmission between the two leads. 
From the mathematical point of view, the problem of 
non-equilibrium OC can be reduced to a non-abelian 
Riemann-Hilbert problem^!, which has not been solved 
in the general case. Progress has been made in the limit 
of long times, where the OC is characterized by the com- 
bination of a power-law d ecay with new exponents and 
weak exponential dampin g ^ 20 ! 21 !. 

We now argue that it is possible to study non- 
equilibrium OC in cold atomic gases by generalizing the 
setup proposed above to the case of a multi-component 
Fermi gas. For simplicity, let us consider a two- 
component gas, with pseudospin species \u) and \d). Our 
goal is to realize a situation in which two components are 
at different chemical potentials (playing the role of two 
leads), and an impurity that can scatter the fermions be- 
tween two species. The first condition can be achieved by 
preparing an imbalanced Fermi gas, with different Fermi 
energies of two components, E F = Eo, E F = Eq + A/i. 
The second condition is more difficult to attain: an im- 
purity's scattering matrix is diagonal in |u), \d) basis, 
with phase shifts at the Fermi energy given by S u _d- In 
the solid state analogy, this corresponds to an impurity 
always fully reflecting electrons; effectively, this brings us 
back to two copies of the equilibrium OC. 

To overcome this difficulty, we consider an application 
of 7r/2 pulse to the pseudo spin of the host fermions. 
After that, the |1) = ^B- and \2) = states 
of fermions will be at different chemical potentials. Cru- 
cially, the impurity's scattering matrix is non-diagonal in 
|1), 1 2) basis (it is obtained by a rotation of the S-matrix 
in the \u), \d) basis by matrix 1+ J^ y )■ The |1), \2) species 
play the role of the electrons in the left and right leads 
in the mesoscopic experiment, with the impurity being 
characterized by a non-trivial scattering matrix. This is 
exactly the situation needed for realizing non-equilibrium 
OC. Performing the Ramsey or spin-echo experiments on 
the impurity pseudospin then allows one to study the re- 
sponse of the non-equilibrium Fermi gas. 

Looking beyond OC, the analogy between the setup 
we just considered and a quantum point contact (QPC) 
with an impurity controlling the transmission through 
the QPC suggests the possibility of studying the full 
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counting statistics of charge transfer. Although theo- 
retically charge counting statistics has played an impor- 
tant role^l, its experimental studies in mesoscopic sys- 
tems have been quite limited^. Simulating quantum 
transport in cold atom experiments is also attractive 
because time-of-flight experiments allow energy-resolved 
measurements; this can be used to study correlations be- 
tween number of transmitted particles at different ener- 
gies. 



VI. CONCLUSIONS 

In conclusion, we studied universal OC in cold atomic 
systems, and discussed the manifestations of OC in the 
Ramsey and spin-echo interfcrometry, RF spectroscopy, 
as well as in energy counting statistics. Beyond the 
equilibrium OC, we have proposed a set of experiments 
which probe the dynamics and transport in nonequilib- 
rium Fermi gases. This provides a useful connection be- 
tween the cold atom physics and mesoscopic physics. 

The basic ingredients required for the experimen- 
tal implementation of our proposal are the following: 
(a) a quantum degenerate Bose-FermP^2U or Fermi- 
FermPSSH j^j^yj-g^ ^he ability to trap one type of 
atom by a strong optical lattice potential, (c) the control 
of the interaction strength between the impurity and host 
atom via changing the impurity hyperfine state, and (d) 
the ability to achieve temperatures that are low enough 
to observe the OC. 

We emphasize that all four requirements are achiev- 
able with currently available techniques, thus we expect 
that our proposal can be implemented in the near future. 
Below we describe some relevant experiments, which, we 
hope, will help to identify the systems which are most 
suitable for studying OC. 

(a) Various quantum degenerate mixtures have been 
realized by multiple group^SHUl 

(b) Recently, the localization of minority atoms by an 
optical lattice was demonstrated with an unbal- 
anced Bose-Bose mixture of 87 Rb and 41 K atoms^. 
The localization of the impurity atoms at length 
scales of roughly 10% of the Fermi wavelength has 
been achieved at typical densities; thus, the im- 
purity can be treated as point-like, and one can 
neglect its excitations to the excited states of the 
trapping potentials, as in our analysis above. 

(c) RF pulses have been used to switch between the hy- 
perfine states of the impurity atom, which interact 
differently with the host fermionsSI. This should 
enable the Ramsey and spin-echo type experi- 
ments that reveal OC. Furthermore, experiments 
with strongly imbalanced Fermi-Fermi mixtures, 
which addressed the polaron dynamics (mobile im- 
purities }p3iH2] demonstrated that the impurity-host 



atom interactions can be tuned in the full interac- 
tion range, from strong attractive to strong repul- 
sive regime. This should allow the exploration of 
different regimes of the OC, discussed above. 

(d) In current experiments, temperatures as low as 
a few percent of the Fermi temperature can be 
achieved 46 * 69 * 85 !. As follows from our analysis 
above, this should be sufficient to observe both the 
universal power law decay as well as the "bottom 
of the band" oscillations with period 2%%/Ef- In 
RF absorption spectra these oscillations result in a 
cusp- like singularity at the Fermi energy Ep. 

It is also worth noting that, for very short times our 
results may be relevant to very heavy mobile impurities 
(mass much larger than that of the host atoms). OC 
will manifest itself, e.g., in the RF spectroscopy exper- 
iments, which are commonly used to probe mixtures of 
cold atoms. A suitable mix tures with large mass ratio 
are 40 K/ 41 K immersed in ^^Mimmn and i73yb/ 174 Yb 
immersed in 6 Li 4 ^'. 

An alternative experimental route to accessing the 
non-equilibrium response of the Fermi gas to a suddenly 
introduced impurity is to create a local scattering po- 
tential by a narrow laser be am. This has been demon- 
strated in recent experiment d 86 * 87 l While in such a setup 
one cannot perform the RF absorption or Ramsey in- 
terference experiments, an observation of OC should be 
possible through energy counting statistics, see Sec. |IV| 
Furthermore, time-of-flight experiments would reveal the 
full distribution functions of scattering processes, provid- 
ing new information about the non-equilibrium state of 
fermions. 

Finally, we note that the ideas presented above are not 
limited to the case of free fermions; one particularly in- 
teresting extension concerns the case of ID interacting 
fermions in an optical lattice. OC in Luttinger liquids 
(LLs) was predicted to exhibit strong deviations from the 
non-interacting case, showing new power laws depend- 
ing on the strength of interactions, universal power-law 
asymptotic behavior, as well as new scaling laws at inter- 
mediate times®'. To access the possibility of observing 
these phenomena in a ID optical lattices, we have carried 
out a numerical studjJ^ of interacting spinless fermions 
in ID lattice; we found that the modifications of OC in 
LL can be observed in finite optical lattices with only 
^100 atoms, which put them within the reach of cur- 
rent experiments. 
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Appendix A: Details on the numerical procedure 

In order to obtain the universal overlap function S(t), 
we evaluate determinants of type (11) numerically. To 



this end, we consider a finite system confined in a sphere 
of radius R whose eigenstates are in the absence and in 
the presence of the impurity 



4>n{r) = \ n sin(fc„r) , k n R = nn 



spectra A(uj) exhibits power law decays and band edges, 
which renders a numerical evaluation of the Fourier trans- 
form extremely challenging, since S(t) has to be known 
for extremely long times. To circumvent this problem, 
we calculate S(t) exactly from Eq. (Ill up to a certain 
time t < t* ~ bOOHEp 1 and then continue S(t) to longer- 
times by fitting it to its asymptotic form (12 1 and (151, 



respectively. Even though, with that we have access to 
S(t) for long times, the Fourier transform is still finite 
dimensional, which gives rise to wild oscillations at the 
edges of A{uj), known as Gibbs phenomenon. We reduce 
the Gibbs phenomenon by applying a Lanczos filter to the 
Fourier coefficients. With this procedure, we produce the 
universal RF-spectra shown in Fig.|4j Importantly, the 
positions of the band edges in the spectra match exactly 
the outcome of Eqs. (fl6j) and (fl7|). 



Appendix B: Universal radio- frequency spectra 



= A n \ — sin(k' n r + 5 n ) , k' n R + S„ = irn 



(Al) 



respectively, where A n = Vyl + ^fnf an d S n — 

— tan _1 (fc^a). For a > there is also a bound state which 
must be treated separately. In principle the evolution of 
the determinants requires even at zero temperature mul- 
tiplications of infinite dimensional matrices. However, 
truncating the infinite number of intermediate states still 
allows to evaluate the determinants with very high accu- 
racy, since the OC is determined by low energy physics. 
In order to obtain the universal overlap functions S(t) 
from the finite system we take a finite size scaling by 
keeping the density constant while taking the system size 
to infinity. Importantly, in a finite size system the overlap 
function S(t) does not decay to zero but rather exhibits 
revivals after sufficiently long times characterized by the 
Fermi velocity and the system size — an additional aspect 
that should be observable in experiments. 

We evaluate the universal RF-spectra A{uS) from the 
Fourier transform of the overlap function S(t). The RF- 



Here we provide the analytic formulas that describe the 
new feature in the absorption spectra A(u>) at frequency 
%uj\ = AE + E F for a < and huj 1 = AE + E b for a > 0, 
the origin of which is illustrated in Fig. 2jjc) (a < 0) and 
(f) (a > 0). The configurations which lead to this non- 
analyticity are given by the excitation of a fermion from 
the bottom of the band to the Fermi surface. The nu- 
merical results in Fig. [4] indicate that the non-analytic 
feature develops into a stronger peak as unitarity is ap- 
proached. This phenomenon, as well as the full structure 
of A(uj) at kp\a\ 3> 1 can be understood as a simple 
interplay between two-body physics happening near the 
bottom of the band, and physics of multiple particle-hole 
excitations being created at the Fermi surface. 

In the time domain, the contributions from the Fermi 
surface excitations and the dynamics of the hole at the 
bottom of the band factorize. The former contributions 
manifest as the usual power laws with unitary phase shift, 
while the latter can be simply evaluated within two-body 
theory, see^3 for details. In the frequency domain, the 
result, written as a convolution of the two terms corre- 
sponding to these two processes, reads: 



tvi T^TB /*mm / HZ° [U-U1 + - )(w-wi + -) F(e) + reg., (Bl) 



T{\/A){E F /hy/± 2tt V 1 hJ \ 1 h 
where F(e) is the probability to excite a hole with energy e measured from the bottom of the band, 

tt-^t + e k F a 

2ma A 



Numerical prefactor is evaluated similar to C, and func- 
tion F(e) accounts for the existence of the bound state 
on the repulsive side. Terms denoted as reg. account for 
a regular contribution at w\ which is not singular at large 



k F \a\. 

Although for < e <ti h 2 /2ma 2 function F(s) be- 
haves as oc y/e, and leads to a weak non-analyticity, for 
h 2 /2ma 2 <C e, it behaves as oc I/a/s- Right at unitarity 
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the scale disappears, and one obtains in A(ui) the di- 
vergence with the universal exponent 1 /4 and a universal 
shoulder ratio, see Eq.(18), which is valid for \k F a\ = oo 
and h\uj — uji \ <C Ep. For large but hnite negative k F a, 
this peak gets smeared out at energies of the order of 
Ti 2 /2ma 2 from its maximal value, as is illustrated in the 



inset to Fig. |4|a). On the repulsive side, the true bound 
state with energy E\> "pinches off' from the bottom of 
the band and leads to a threshold with an exponent 3/4. 
Thus for large but finite kpa > the universal form of 
A(u>) near lo\ has a characteristic double peak structure, 
as is seen in Fig. gb) for k F a = 6.0. 
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